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Abstract 



O 

o 

P^ . We consider network models of quantum localisation in which a particle with 

a two-component wave function propagates through the nodes and along the 

edges of an arbitrary directed graph, subject to a random SU(2) rotation on 
CS| ' each edge it traverses. The propagation through each node is specified by 

an arbitrary but fixed 5-matrix. Such networks model localisation problems 

^ . in class C of the classification of Altland and Zirnbauer, and, on suitable 

^T , graphs, they model the spin quantum Hall transition. We extend the analy- 

(^ \ ses of Gruzberg, Ludwig and Read and of Beamond, Cardy and Chalker to 

^O ' show that, on an arbitrary graph, the mean density of states and the mean 

conductance may be calculated in terms of observables of a classical history- 
^D ■ dependent random walk on the same graph. The transition weights for this 

,^ , process are explicitly related to the elements of the 5-matrices. They are 

correctly normalised but, on graphs with nodes of degree greater than 4, not 
T^ ' necessarily non-negative (and therefore interpretable as probabilities) unless 

C^ ' a sufficient number of them happen to vanish. Our methods use a supersym- 

metric path integral formulation of the problem which is completely finite and 

rigorous. 

•i-H 

X 

C3 ■ I. INTRODUCTION 

Network models of quantum localisation were first introduced by Chalker and Coddington 
[1] to model the transition between plateaux in integer quantum Hall systems. Reduced to 
their essentials, they describe the propagation of a single quantum-mechanical particle along 
the directed edges and through the nodes of a graph. For the Chalker-Coddington model, 
this graph is some large but bounded domain of the L-lattice, a square lattice whose edges 
are directed in such a way that the particle turns through ±90° at each node. In propagating 
along each edge, the single-component wave function is multiplied by random phases, which 
are i.i.d. random variables with a uniform distribution in [0,27r). The propagation through 
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each node is described by a unitary 2x2 S'-matrix of amplitudes between the two incoming 
and two outgoing edges. 

The integer quantum Hall plateau transition is only one among several possible uni- 
versality classes of quantum localisation transitions, which have been classified by Altland 
and Zirnbauer [2] according to the symmetry properties of the underlying single-particle 
hamiltonian Ti. Another, known as class C, corresponds to the existence of a symmetry 

(TyT-C(jy = — 7i* . (1-1) 

This gives rise to a pairing between eigenstates with energies ±E, while E = is special 
and may correspond to delocalised eigenstates, even in two dimensions. Class C is supposed 
to be realised in disordered spin-singlet superconductors in which time-reversal symmetry is 
broken, but Zeeman splitting is negligible [2]. The fact that spin is still conserved can then 
lead to a spin quantum Hall effect. 

The appropriate network model for this on the L-lattice was formulated and studied nu- 
merically by Kagalovsky et al [3]. An equivalent spin-chain hamiltonian was also investigated 
by Senthil et al [4] . However, in a remarkable paper, Gruzberg, Ludwig and Read [5] showed 
that the mean single-particle Green function, as well as the mean conductance, may be ex- 
pressed in terms of classical averages of appropriate observables of the hulls (boundaries) of 
clusters in classical bond percolation on the square lattice. The critical exponents [5] and 
some other universal properties [6] of the spin quantum Hall transition in two dimensions 
are thus exactly known. 

The methods of Gruzberg et al [5] used supersymmetry to average over quenched disorder. 
One of the essential features of this, which will also appear in our analysis, is the reduction 
of the Hilbert space on each edge to one of finite dimension. They then analysed the transfer 
matrix for the L-lattice, and demonstrated its equivalence to that for percolation hulls. 

One of the interesting features of percolation hulls on the square lattice is that they may 
be generated, independently of the underlying percolation problem, as history-dependent 
random walks on the L-lattice. Consider a random walk which begins on some edge, and 
steps through one node in unit time. At each node it turns to the left or right with prob- 
abilities p or 1 — p respectively. However, it cannot traverse a given edge more than once, 
so that whenever it returns to a node it has already visited it is forced to exit along the 
other (empty) edge. Eventually, on a closed graph, it will return to its initial edge. The 
statistical properties of such loops are identical to those of a single closed hull in the per- 
colation problem. On the L-lattice in the thermodynamic limit, when J9 7^ |, loops close 
almost surely after a finite number of steps. This corresponds to quantum localisation in 
the network model. The only delocalised states occur when p = |, at the bond percolation 
threshold. 

Motivated by this work, Beamond, Cardy and Chalker [7] investigated class C network 
models on arbitrary graphs. Their methods did not use supersymmetry. Instead they showed 
that, for quantities like the average Green function G and the mean conductance (related 
to |Gp), there is a massive cancellation between paths in the Feynman expansion of these 
quantities which leaves essentially classical paths whose weights correspond to those of a 
history- dependent random walk. This reproduces the results of Gruzberg et al [5] when 
specialised to the L-lattice. 



However, the proof of Beamond et al [7] was restricted to graphs in which each node has 
N = 2 incoming and outgoing directed edges. It becomes too cumbersome to generahse it to 
graphs with nodes with N > 2. Nevertheless, it seems important to find such a generahsation 
in order to be able to investigate, for example, the properties of such network models on 
simple regular lattices embedded in three or more dimensions. 

In this paper, we find this generalisation. Like Gruzberg et al [5], we use supersymmetry 
to perform the quenched average. However, we do this within a path integral, rather than 
a Hilbert space, formulation of the model, which allows for the treatment of an arbitrary 
graph, not only those regular lattices which admit a transfer matrix. The result is positive, 
in the sense that we can prove that the mean of G and of IGp may be expressed as a sum 
of history-dependent classical random walks on the graph. The weights at each node, given 
that a given set of incoming edges is occupied by the walk, when summed over all possible 
outcomes, correctly sum to unity. However, in general these weights are not positive, and 
therefore cannot be interpreted as probabilities. In fact we show that, for N > 2, this 
condition can only be satisfied if a certain number of elements of S vanish. A sufficient 
condition for this is that the ^-matrix at each node is a direct product of S'-matrices for 
2 — !► 2 nodes, that is, the node may be decomposed into 2-^2 nodes (in which case, of 
course, the analysis of Beamond et al [7] applies.) For A^ = 3 and 4 we have shown that this 
condition is also necessary. 

Despite this negative result, the proof of the general theorem sheds further light on the 
analysis of Gruzberg et al [5], as well as giving a more elegant derivation of Beamond et 
al [7]. These methods may also be used to determine which combinations of averages of 
higher-point Green functions in the network model may be related to observables in the 
classical problem, thus giving a simpler derivation of some of the results of Mirlin, Evers, 
and Mildenberger [8]. 

The layout of this paper is as follows. In the next Section, we define the network models 
and the observables of interest. Then we are able to state our main Theorems (1 and 2). In 
Sec. HI we introduce the path integral machinery necessary to compute these observables, 
and to perform quenched averages. The supersymmetric path integral involves both bosonic 
(commuting) and fermionic (anticommuting) variables on each edge, but after the quenched 
average is taken we show that these reduce to the propagation of only single fermion-fermion 
or boson-fermion pairs. In Sec. IV we then show that the propagation of a fermion-fermion 
pair through the lattice obeys the rules of a classical history-dependent random walk, once 
all the other degrees of freedom are traced out. This proves Theorem 1. We also consider the 
case of open systems, and conductance measurements. Finally, in Sec. V, we consider the 
probabilistic interpretation in terms of history-dependent random walks, and prove Theorem 
4 which states conditions under which these weights are non-negative. 

II. DEFINITION OF THE MODEL AND OBSERVABLES. 

Let ^ be a graph consisting of J\f directed edges, and nodes. Initially we consider only 
closed graphs. At each node there is an equal number of incoming and outgoing directed 
edges. Apart from this, Q is arbitrary. We wish to define a network model on this graph with 
describes the propagation of a quantum-mechanical particle whose hamiltonian 7i obeys 



the symmetry (1.1). First we note that the single-particle Hilbert space must be even- 
dimensional in order to be able to define the action of (jy. For simplicity we take this to be 
two-dimensional (in [7] a method was proposed for reducing any class C network model with 
an even-dimensional single-particle Hilbert space to this case.) Since the network describes 
propagation over finite time steps At, we need to define the unitary evolution operators Ue 
and Un which evolve the wave function along each edge and node respectively. By (1.1), 
they must obey CyUcTy = U* . Each Ue, therefore, must be an element of Sp(2)~SU(2). We 
take these to be i.i.d. random variables, each uniformly distributed with respect to the Haar 
measure of SU(2). Quenched averages with respect to this measure will be denoted by an 
overline. f/„ is the unitary S'-matrix for the node n. As will become clear, it suffices to take 
this to be diagonal in the SU(2) indices, so, by the above, it must be real and therefore an 
element of 0(iV), where A^ is the number of incoming (and outgoing) edges at the node. 
The Un may vary from node to node, and in principle they may also be random variables. 
However, our theorems apply to a fixed realisation of these S'-matrix elements. 

The full unitary evolution operator U for the whole network is then a direct sum over 
the edges and nodes of 1 ... f/g ... ® 1 and 1 (g) ... f/„ ... eg) 1. Of interest is the 
Green function, which is the matrix element of the resolvent operator (1 — zU)~^ between 
states localised on two edges: 

G(e2,ei;z) = (621(1 -^W)-^|ei) . (2.1) 

This is a 2 X 2 matrix in SU(2) space. For |z| < 1 this may be expanded as a sum of Feynman 
paths on Q beginning at Ci and ending at 62: each path is weighted by an ordered product of 
SU(2) matrices Ue along the edges is traverses, and a product of S'-matrix elements according 
to how it passes through each node, as well as a factor z raised to the power of its length. 
However, each edge may be traversed an arbitrary number of times. If we were to formulate 
the sum over such paths using a transfer matrix formalism (assuming that Q allows this), 
the Hilbert space on each edge would be infinite-dimensional. 

Alternatively, for \z\ > 1, we may write the resolvent as —z^^U'' {1 — z^^W)^^ and expand 
G in powers of z^'^, as a similar sum over paths. Each path is now weighted by an ordered 
product of factors z^^U^, as well as an overall factor of —z~'^. 

The Green function may be used to compute the density of states of U via its diagonal 
elements G{e, e; z). For a closed graph, these are of the form exp(iej). We define the density 
of states p(e) = J2j ^(e ^ ^j)- Then 

Pie) = ^^ E li^ (Tr G{e, e; (1 - 6)e~n " Tr G{e, e; (1 + 6)6'^) (2.2) 

An open system may be defined by cutting open a subset of the edges of Q. We may then 
perform a conductance measurement by attaching leads to a subset {cin} of the incoming 
edges, and to a subset {cout} of the outgoing edges. The transmission matrix t between 
these two leads has elements (eout|(l — W)|ein), and the conductance is then given by the 
Landauer formula as g = Tr tH. In particular, the point conductance between two edges is 
TrG'(eout,ei„; l)^"G'(eout, e^; 1). 

The above defines the quantum problem which we wish to study. Our main theorems will 
relate it to a classical problem, defined as follows. For each node in Q, adopt an arbitrary 



but fixed labelling of the incoming edges j G {!,..., A^} and outgoing edges i E {1, . . . , N}, 
and denote the elements of the corresponding S*- matrix by Sij. Note that det S* = ±1, with 
the sign being dependent on the choice of labelling. 

Define a trail r on ^ as a sequence of distinct edges (ei, . . . , e\r\) such that e^ and e^+i 
are incoming and outgoing edges of the same node, for each 1 < k < \t\ — 1. It is a (rooted) 
closed trail if e\r\ and ei also share the same node. Note that a trail cannot pass along a 
given edge more than once, but it may pass through a given node any number of times, up 
to its order. 

For a particular trail r, and a particular node n, denote the set of incoming edges on r by 
Jn;T, and the set of outgoing edges by /„.^. These are both ordered subsets of {1,2,..., N}. 
A given trail associates an element of In-r with each element of Jn;T, and vice versa, and 
thus may be associated with a permutation iTn^r of the ordered elements of Jn-r- Denote the 
signature of this by (—1)'^"'^. Let det Sj^j denote the minor of S restricted to the ordered 
subsets I and J of the outgoing and incoming channels respectively. 

We are now ready to state 
Theorem 1. The mean of G{ei,e2; z) vanishes if ci ^ 62, while in the case of equality it is 
given by 

f 2 — y^ . X 7;W ^ z2|r(e)| . |r| < 1 

TrG'(e,e;^) = ^ ^-(^)'^^!^2 4)| ' ! (2-3) 

where the sums are over all closed trails r(e) rooted at e and Wr(e) is the weight of each, 
given by the product over all the nodes on r(e) of factors 

fi(/„^,;J„^,)^(-l)-- n ^-.;.0-M (det 5,„^.,j„ J (2.4) 

Note that the first two factors are the term in the expansion of det Si^.^^j^.^ corresponding 
to the permutation 7r„;r: if we were to sum (2.4) over all permutations, we would obtain 
(det S'/„.^^j„.^) . (2.4) is a generahsation of the main result of [7], which applies to the 
case N = 2. In this case, the elements of S may be taken to be 6*11 = S22 = cos9 and 
S^2 = -S21 = sin^. If J = {1} and / = {!}, or if J = {2} and / = {2}, the weight is 
cos^ 9. If J = {1} and / = {2}, or if J = {2} and / = {!}, it is sin^ 9. But if J = {1, 2} and 
/ = {1,2} or {2,1}, it is unity. 

The next theorem gives the equivalent result for conductance measurements. 
Theorem 2. The mean point conductance g between two edges Cin and Cout is given by twice 
the sum over all open trails on Q connecting the two edges, each such path being weighted as 
for the closed loops in Theorem 1. 

Thus, if the weights on the trails can be interpreted as a probability measure, the mean 
conductance between two contacts is just twice the expected number of open trails which 
connect them. 

The weight for a single trail r is given by a product of weights corresponding to each 
node thorough which r passes, once the whole of r is given. Alternatively, we may build 
up these weights as a product of factors incurred each time r passes through a given node. 
For example, the first time it passes through, entering via edge ji and leaving by edge ii it 
incurs a weight, according to (2.4), of VL{ii]ji) = Sf_^j_^. If it passes though the same node 



again, this time entering along j2 and leaving along i2, it then incurs a conditional weight^ 
U!{ii,i2;ji,J2) = ^{ii,i2;ji,J2)/^{ii;ji) and so on. In general 

wit,...,t,;n,...,j^ = ^^'^-r''-'''''^'---^''''-^^ (2.5) 

i2(?i...,2p_i;ji,...,jp_i) 

The next theorem states that these conditional weights are properly normalised, in the sense 
that they give unity when summed over all possible outcomes: 
Theorem 3. The weights w{ii . . . ,ip;ji, . . . , jp) satisfy 

^ w{ii...,ip;ji,...,jp) = 1 (2.6) 

ip({ii,...,ip-l} 

Thus, as long as they are non-negative, they define a set of transition probabilities for a 
discrete random process whereby the ensemble of trails may be dynamically generated with 
the correct weights, the trail growing by one unit at each time step. Since the weights at a 
given node depend on whether (and how) it has been visited in the past, the process may 
be thought of as a history-dependent random walk. When all nodes have N = 2, this is 
straightforward: the first time r passes through a given node, it incurs a factor cos^ 6 or 
sin^ 9. If it passes through a second time, this factor is unmodified. 

We may ask whether this positivity can extend to nodes with A^ > 3. To answer this, 
we need a notion of reducibility. The S'-matrix at a node with A^ > 3 is said to be reducible 
if it admits a factorisation of the form S = S^^^S^'^\ where (after a possible re-ordering of 
the incoming and outgoing channel labels) the N x N S^^^ and S^"^^ matrices have the block 
diagonal forms 

where s^^ and s}^_ are orthogonal p x p and (A^ — Q') x {N — q) matrices respectively, and 
p > q. This is illustrated in Fig. 1. This procedure may be repeated. An N x N S'-matrix 
is said to be completely reducible if it can be factorised in this way into 2^2 S'-matrices. 
Theorem 4. At a node with N > 3, a sufficient condition for the weights (2.4) to be all 
non-negative is that the S-matrix is completely reducible. For N = 3 and 4 this is also 
necessary. Thus, in these cases, the network model on Q could have been described on an 
equivalent graph with only N = 2 nodes. 

III. PATH INTEGRAL REPRESENTATION. 

In the standard way, G may be written as a path integral over commuting (bosonic) 
variables. The notation is a little complicated, but the basic idea is simple. Label each end 



^This assumes 5'jijj 7^ 0. If the conditional weight can be interpreted as a probability, 
'w{ii,i2',ji,J2) ^ 1 and therefore it has a finite limit as Si^^j-^ -^ 0. Even if this is not the case, the 
unconditional weight w{ii,i2; ji, J2) wiii', ji) vanishes in this limit. 




FIG. 1. A reducible S'-matrix. 

of a given directed edge e by e^j and e^, in the direction of propagation e^i ^ e^. Introduce 
complex integration variables 6i?(e) and biie), which are each 2-component column vectors in 
SU(2) space, their components being labelled bjiaie) and hLa{,&) respectively, where a = 1,2. 
Then G can be written as 



^(62,61;^) = (&L(e2)6l(ei)) 
where Wh = Wedee + W^node with 



JUe[dbLmdbRie)]bLie2)biie^)e'^>' 



lUe[dbL{e)][dbR{e)]e 



w. 



W,dse = zJ2b{{e)UMe) 

e 
W^ode = Y.HY.^*Ra{^i)^Sn)ijbLa{,ej] 



n a ij 

and the integration is wrt the usual coherent state measure 



f[db] = (l/vr^) fe-^^^l[dRebadlmba 



(3.1) 

(3.2) 
(3.3) 

(3.4) 



Note that there is a finite number of integrations, if Q is finite, and that no time-ordering 
necessary: we can imagine writing everything out in terms of components, and all quantities 
in the path integral are commuting. On a finite graph, only a finite number of integrations 
is necessary. The exponentiation of Wb correctly takes into account the multiple traversing 
of edges by Feynman paths. 

The next step is to average over the quenched random variables Ue- As usual, since these 
occur in both the numerator and denominator, this is most easily done either by introducing 
replicas, or by adding an anticommuting (fermionic) copy of the bosonic variables, making 
it super symmetric. We opt for the latter. Thus to each pair of complex integration variables 
{b\ b) we introduce a pair of Grassmann variables (/, /) with corresponding labels, and we 
add to Wb a term Wf of identical form with bosonic variables replaced by fermionic ones. 
The Grassmann integration is defined by 



[df] 



dfdfe 



Jf 



(3.5) 



so that 



J[df]f = J[df]f = 0; (3.6) 

J[df]l = J[df]ff = l (3.7) 

Integrating over the fermionic variables cancels the denominator in (3.1), so that we may 
write 

G{e2,e^) = /nM&L(e)]M6ij(e)]M/L(e)]M/^(e)]6^(e2)6l(ei)e^^+^/ (3.8) 

•^ e 

However, Wb + Wf is invariant under global supersymmetry, so G may equally well be 
expressed, for example, as (/L(e2)/L(ei)). 

Quenched average. 

The average over the SU(2) matrix f/ on a given edge has the form 

J dU exp{zb{Ubn + zhUfR) (3.9) 

where the integral is with respect to the invariant measure on SU(2), normalised so that 

JdU = 1. 

Lemma 1: The above integral equals 1 + 1^;^ det M, where M is the 2x2 matrix with 

components Mij = bl^Rj + fufRj- 

Proof: Because the fermionic variables have only two components, and any such component 

squares to zero, the expansion in the fermionic part terminates: 

JdU exp{zb[Ubn) (l + zhUfn + \z\hUfRf) (3.10) 

The first term, the purely bosonic integral, is identically equal to unity. This follows from 
the observation that the integral is invariant under the substitutions bji -^ XVubji, b\ — >■ 
\~^b\yl^ where Vl and Vr are independent SU(2) matrices, and A is a complex number, 
and there is no combination of b\ and bji which has this property, save a constant. However, 
an explicit proof is given in Appendix A. 

The third, purely fermionic, term is also easy: 

{hUfR? = ihlUnfm + /L2f/2l/m + /Llf/l2/i?2 + fL2U22fR2? (3.11) 

= 2{fLlUufRl){fL2U22fR2) + 2{fL2U2lfRl){hlUl2fR2) (3.12) 

= 2hlfL2fR2fRl{UllU22 - U,2U2l) (3.13) 

= 2hjL_2fR2fRl _ (3.14) 

= detf(^^{«^ ^^'i""') (3.15) 

V JL2JR1 jL2jR2 / 

where the fourth line follows because det U = 1. The expression is therefore independent of 
U, and the integration is then the same as in the purely bosonic term, which gives a factor 
1 as before. 



The second term can also be worked out explicitly, but it easier to invoke the supersym- 
metry, and simply add ^^j^rj to each element fnfRj of the above matrix. Note that the 
purely bosonic part of the determinant vanishes. I 

The result of the quenched average over the SU(2) matrix on a given edge is therefore 

1 + ¥\blJL2 - bljLi)ibnifR2 - bR2fm) + z'ihjL2)ifR2fm) (3.16) 

The interpretation of this is clear: after averaging over the SU(2) matrices, the only paths 
which contribute are those in which on each edge the allowed propagation is of either the 
identity, a pair of fermions /1/2, or a boson-fermion pair (l/v/2)(&i/2 — /i&2)- Note that in 
each case the combinations in parentheses above are SU(2) singlets. Note also that, having 
averaged over the edge variables Ue, the distinction between L and R is now immaterial, 
and we can henceforth drop these labels. 

The above result has several important consequences. First, there is now only a finite 
number 3 of possibilities for propagation along the J\f edges of a finite graph Q. (This is 
equivalent to the result of Gruzberg et al [5] that the Hilbert space of the transfer matrix 
is finite-dimensional.) Second, it is clear why the assumption that the scattering at the 
nodes in diagonal in the SU(2) indices was not crucial: only the singlet invariant amplitude 
matters. Third, the only non-zero two-point functions with 62 7^ ei are 

i((&i(e2)/2(e2) - &2(e2)/i(e2)) {bl{e,)Me^) - b;ie,)J,{e,))) (3.17) 

= (/i(e2)/2(e2)72(e i) A(ei)) (3.18) 

= G'iiG'22 - ^12^21 = det G{e2, ei; z) (3.19) 

Let us for the moment take z real. Then G{e2,ei; z), as a sum over Feynman paths, 
is a linear combination of SU(2) matrices with real coefficients. Any such 2x2 matrix is 
itself proportional to an SU(2) matrix, up to a real scalar (see Appendix.) Thus we may 
write G = \G where A is real and G G SU(2). Hence det G = \^, and G^G = A^J, so that 
Tr G'^G = 2 det G. The right hand side is a polynomial in z. For general complex z we have, 
therefore. 



2 det G{e2, ei; z) = Tr G{e2, ei, z*yG{e2, ei, z) (3.20) 

When z = 1 this is the mean point conductance, which is therefore given, up to a factor 2, 
by the two-point functions in (3.17). 

Since only SU(2) singlets now propagate, it follows that two-point functions like 
(/a(e2)/a(ei)) = G{e2,ei) vanish if 62 7^ ei. This is because, once the matrices Ue have 
been traced out, the supersymmetric path integral possesses a local SU(2) gauge invariance 
under (Ke),/(e)) ^ (V;6(e), V;/(e)) with V; G SU(2). 

However, this does not apply if 62 = ei. In fact, because of (3.17), it follows that 



G(e, e; z)n = (/i(e)A(e)) = {h{e)f2{e)f^{e)f,{e)) = det G(e, e; z) (3.21) 



IV. PROPAGATION THROUGH THE NODES. 

In the last section, we showed that, for a graph Q with J\f edges, the quenched average 
of the path integral can be written as a sum of 3 terms, according to which of the three 
terms in (3.16) (corresponding to the propagation of a bf pair, an // pair, or the identity) is 
chosen on each edge. Let us now consider just one of these terms, and one particular node. 

The contribution to the path integral from this node has the form 

llAaXn)Sll^.^irj) (4.1) 

i 3 

where A^ = 1, A2 = /1/2 and A3 = (1/^2) (61/2 - 62/1) and 

S = exp(^ J2ib:aS^,b,a + LS.jfja)) (4.2) 

a=l i,j 

In doing this, we have brought together in the path integral all the factors associated 
with the given node. There is a subtlety, however, because the boson-fermion variables A3 
and A3 anticommute with each other. At a given node, we may arrange these factors in 
the standard order determined by the fixed (but arbitrary) labelling of the incoming and 
outgoing edges. For a given term out of the 3 possibilities, this will introduce an overall 
factor ±1. 

Define a decomposition of the node as a pairing of each outgoing edges i with a unique 
incoming edge j. This defines a permutation n of the edge labels, whereby the outgoing 
edge i paired with the incoming edge j is vr(j). Carried through for every node in turn, this 
decomposes Q into a union of disjoint directed closed loops (and open paths if Q is open), 
such that every edge lies on just one loop or open path, and each loop or open path may pass 
along a given edge no more than once. The following Proposition shows that we are allowed 
to do this inside the path integral, as long as we weight each decomposition correctly: 
Proposition 1. The result of performing the integration over the variables {bj,fj) and 
{blifi) in (4-1) is the same as if S were replaced by 

det S E(-l)" I[knij)S^,S^^,^^ (4.3) 

TT ij 

that is, it is given by a weighted sum over all decompositions tt. In each decomposition, each 
state on the incoming edge j propagates freely to vr(j). 

Proof: Since the numbers of each component of both bosons and fermions are the same in 
the incoming and outgoing channels, and bosons are always paired with fermions, it follows 
that the numbers of // and fb pairs are individually conserved at every node. Let us call the 
subsets of the A^ outgoing channels occupied by an // pair, a bf pair, or empty, FF, FB, and 
E, respectively, and similarly for the incoming channels, FF, FB and E. The integrations 
may now be performed, expanding S to second order in the Sij and using Wick's theorem. 
Each fermion (boson) in outgoing channel i, when contracted with a fermion (boson) in 
the incoming channel j, gives (up to a sign) a factor SabSij. The bosons in FB may only 
contract onto the bosons in FB, but the complication is that some of the fermions in FF 

10 



may contract onto fermions in FB, and some of those in FB may contract onto those in FF. 
However, every set of possible contractions will involve each outgoing channel in FF U FB 
and each incoming channel in FF U FB exactly twice. Thus, if a denotes a permutation of 
the channels in FF U FB, then, the general form of the result will be 

^aa,a' n ^iAi) n Sii^„i(ii) (4.4) 

a,a' ieFFUFB i'f^FFUFB 

where the a„y are numerical coefficients. 

We have already introduced the notation det Sjj for the minor of S restricted to the 
ordered subsets I and J of outgoing and incoming channels. Now define perm Sij to be the 
corresponding permanent, that is, with all the terms having the same sign +1. Then the 
claim is that the result of the integration is 

ciet >^ pp^^ ' perm o pj^ "pg ■ ciet i^ ppuPB Wfufb \ ^^) 

This expression has the correct properties in that: (a) each channel index appears exactly 
twice in each term; (b) Sij with i G FF and j G FB (and also with i G FB and j G FF 
occurs at most once; (c) it is symmetric under permutations of the channels in FF, and 
separately in FF; (d) it is antisymmetric under permutations of the channels in FB, and 
separately in FB; and (e) it has the correct overall numerical coefficient. 

In order to prove (4.5), it is helpful ffist to consider what happens if each SU(2) singlet 
75(^1/2 — &2/1) is replaced by 61/2 (and similarly for the conjugate variables in the incoming 
channels.) In that case, the result follows immediately. The fu with i G FF can contract 
only onto the /^ with j G FF, giving the ffist factor in (4.5). Similarly, the bu with i G FB 
can contract only onto the blj with j G FB, giving the second factor. Finally the /2J are 
free to contract onto any of the /2J, leading to the last factor. The reason that this result 
continues to hold when each 61/2 is replaced by the singlet combination is the local gauge 
invariance already alluded to: we could imagine multiplying the whole amplitude by an 
independent SU(2) matrix in each channel, and averaging over this. The final result, being 
gauge invariant, would not change, but it would project 61/2 onto the singlet combination. 

We now need the following property of O(A^) matrices: 
Lemma 2. If S E 0{N), and det S" and det S"' are complementary minors of S, then 
det S' = det S" -detS. 

Proof: This relies on the fact that if S' has rank p, and Tj^,,,j^ is a tensor of rank p, then 
^ji--jpjp+i-JN'^ji-jp (where e... is the Levi-Civita symbol) transforms under proper rotations 
as a tensor of rank N — p, and changes sign under parity. I 

In our case this implies that det 5'^i;'uF_B ffufb ~ ^^^ ^e^ ' ^^^ ^1 ^^ ^'^^^ (4-5) reads 

det Sppj^ ■ perm Sp^ -p^ ■ det Sp p ■ det S (4.6) 

Now look at (4.3), inserted into the path integral instead of S. The Kronecker deltas which 
conserve the labels a restrict the sum over permutations n to those which map FF onto some 
permutation npp of FF, FB onto some permutation Tips of FB, and so on. The signature 
(— l)'^ decomposes into a product of the signatures of the three permutations. Now, since 
the // pairs propagate freely (and they commute among themselves), the integrations over 
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these variables give unity. The sum over npF is therefore Y^Trppi—'^Y'''' YiieFF Si^-KpFii) — 
det Sppj^, which gives the first factor in (4.6). Ahhough the bf pairs also propagate freely, 
they are fermionic, which mens that their contractions give rise to an extra factor (—1)'^^^. 
On summing over all irpB, we get the second factor in (4.6). The remaining factors of >S',r(i)j5 
with j G E, when summed over tte, give the last factor. 

We have shown the equivalence of the two expressions S and (4.3) at each node, for each of 
the 3 terms in the expansion of the path integral. We may now restore the anticommuting 
bf factors to their original ordering in the path integral, thus removing the possible overall 
sign. This concludes the proof of Prop. 1. I 

Proof of Theorems 1 and 2. 



First consider the case when Q is closed. In Sec. Ill it was shown that G{e, e) is given 
by the correlation function (/Li(e)/i;,2(e)/L2(e)/Li(e)) in the super symmetric path integral. 
By the results of the previous Section, this is given by a sum of terms in which each edge 
except e is occupied by either an // pair, a bf pair, or the identity (and e is occupied only 
by an // pair.) Moreover, the path integral is given by a sum of terms, each corresponding 
to a decomposition of Q into closed loops. Along all but one of the closed loops can freely 
propagate an // pair, giving an overall factor +1, an bf pair, giving -1, or the identity, 
giving +1. The first two contributions cancel, leaving a factor +1 for each of these closed 
loops. The exception is the unique loop which contains the edge e, which can be thought of 
as a closed trail T(e), rooted at e. Along this only an // pair is allowed to propagate. 

Now sum over all decompositions of Q which contain the specified trail T(e). At a given 
node n, r(e) occupies the incoming edges Jn-r and the outgoing edges In-r- The sum in 
(4.3) includes only those permutations n for which 7i{Jn;T) is some permutation of /n;r- 
This implies that vr acting on the complementary subset Jn-^r is some permutation 7f of the 
complement /„;,-. If we now sum the factors of Sij in (4.3) with i G /„;r and j G J„;t- over 
the permutations tt , weighted by (—!)'", we get det Sj -j . Using Lemma 2 again, this 
equals det Sj^.^^j^.^ ■ det S. The latter factor of det S combines with exphcit one in (4.3) to 
give unity. The remaining factors then give the weight (2.4) of the node n on the trail T(e). 
This proves Theorem 1. I 

Theorem 2 follows similarly. For an open graph, TTG^{e2, ei)G(e2, ei) is given by a sum 
of decompositions of Q as before, into closed loops as well as open paths which connect the 
incoming and outgoing external edges. Along these propagate either // pairs, bf pairs, or 
the identity, with weights at each node given by (4.3). In each decomposition, there is a 
unique open trail r from ei to 62, carrying an // pair. The other open paths must carry 
the identity, otherwise the path integration over the free bosonic and fermionic variables 
at their ends would give zero. They therefore contribute a factor 1. All the other closed 
loops also contribute a factor 1 after the cancellation between the // and bf pairs which 
propagate around each of them. We are left with a single // pair propagating along r. The 
summation over all the decompositions of Q containing a given open trail r then gives a 
factor det Sj^.^j^.^ ■ det S at each node as above. This proves Theorem 2. I 
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V. PROBABILISTIC INTERPRETATION. 

Normalisation. 

We now prove Theorem 3, which states that the weights f2(/, J) in Theorem 1 lead, if non- 
negative, through (2.5) to correctly normalised transition probabilities w{ii, . . . ,ip;ji, . . . , jp) 
for the trail r(e) interpreted as a classical random walk on the edges of Q. A necessary and 
sufficient condition for this is 

J2 ^(^1, • • • 5 ip-i, ip'Ji, ■ ■ -Jp-iJp) = ^{h, • • • , ip-i;ji, ■ ■ • , Jp-i) (5.1) 

ip^{h,...,ip-i} 

Without loss of generality, we may relabel the rows and columns of S so that ik = k for 
1 < k < p — 1, and jk = k ioT 1 < k < p. Notice that we can remove the restriction on the 
sum over ip because the summand formally vanishes whenever 1 < Sp < p—1. The index ip 
occurs on the left hand side of (5.1) in the factor Si ^p as well as in each term of the expansion 
of the minor det "Sii ... jp}.|i ... p}, where it occurs as Si^^k with 1 < k < p. Thus the sum over 
ip, in each term in the expansion of the determinant, has the form J2i Si^^pSi^^k = Spk, 
from the orthonormality of the rows of S. The coefficient of this term is just the subminor 
det >S'{i^...^p„i};{i^...^p_i} which occurs on the right hand side of (5.1). All the remaining factors 
ni<fc<p-i Sk,k are the same on both sides. This demonstrates the validity of (5.1) and thus 
Theorem 3. I 

Positivity of the weights. 

Although we have argued that the weights Q appearing in Theorem 1 are normalised, 
they may only be interpreted as probabilities if they are all non-negative. This places strong 
constraints on the S'-matrix at each node. 

Taking first the case when the sets / and J comprise all the outgoing and incoming edges 
of the node, we see that the weights are all non-negative if and only if every term in the 
expansion of det S has the same sign, or vanishes. In fact, this is also a sufficient condition 
for all the weights to be non-negative when / and J are proper subsets. This is because, 
by Lemma 2, det Si^j is, up to a factor det S* = ±1, the same as its conjugate minor, and 
therefore each term in (2.4) is, up to an overall sign, a sum of a subset of terms in the 
expansion of det S. They therefore all have the same sign, or vanish, if this is true of the 
individual terms in the expansion. 

For N = 2, this is always the case. If det S* = 1, we can write S = [ . „ „ ), 

so that the terms in the expansion are (cos^^, sin^^); or if det S* = — 1 we can write S = 

. „ I, in which case they are (— sin^^, — cos^ ^). However, for an orthogonal 
cos6' -sm6'/ ^ V J ) to 

matrix with A^ > 2, this constraint becomes nontrivial. 
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1. N = 3. 

Consider first the case A^ = 3. It is elementary to show that if the 3! terms in the 
expansion of the determinant of any 3x3 matrix all have the sign (or vanish) then there 
must be at least one vanishing element. For consider the product of all these terms. This 
contains each element Sij exactly twice. There are six terms in all, and three of these, 
corresponding to the odd permutations, occur with minus signs. Hence the product of all 
the terms is — nf=i 11^=1 Sfj < 0. This would be impossible if all the Sij were no n- vanishing. 

Now any 0(3) rotation can be composed of three suitable 0(2) rotations about different 
axes, for example through the Euler angles. This composition may, in general, be pictured 
using a diagram like that in Fig. 2. Each intersection of lines labelled by i and j corresponds 




FIG. 2. Diagram illustrating Euler angle representation of an 0(3) S'-matrix. 

to an 0(2) rotation in the ij plane, represented by an 0(2) matrix s*-"-* with a = 1, 2, 3. The 
element Sij of the full 0(3) matrix is given by a sum over directed paths from j to i in the 
diagram, each path being weighted by a product of the appropriate 0(2) matrix elements. 
For example, 

Si3 = 4? (5.2) 

•Jsi — S32 S21 + S33 S31 Sii [i:).6) 

Each topologically distinct way of drawing and labelling Fig. 2 corresponds to a different 
but equivalent Euler angle parametrisation. 

We can always draw the diagram so that the matrix element which vanishes by the 
above argument (in this example 5*13) is given by a simple form like (5.2). This implies that 
■^13 = -^si = 0) ^^^ therefore that Sii = S33 = 1 (note that s*-^-* can always be chosen as a 
proper rotation.) This means that we can picture the lines 1 and 3 simply crossing at the 
vertex (2), and that the full 0(3) rotation reduces into a product of just two 0(2) rotations, 
as in the definition (2.7) of reducibility. 
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2. TV > 3. 

An 0(iV) matrix has several distinct but equivalent Euler angle representations as a 
composition of ^N{N — 1) 0(2) rotations, which may be pictured using a generalisation of 
Fig. 2. An example for A^ = 4 is shown in Fig. 3. In such a diagram any given line intersects 




FIG. 3. Euler angle representation for A^ = 4: an example where vanishing elements force the 
matrix to be completely reducible. In this case 5i4 = S'13 = 5*24 = 0. 

each of the others exactly once. If the matrix is completely reducible there is at least one 
representation which has a tree structure, that is, contains no cycles. An example is shown 
in Fig. 3. In this case, many elements of S must vanish, and those which do not are each 
given by a single term which is a product of 0(2) matrix elements along a single possible 
path through the diagram. 

Complete reducibility is a sufficient condition for the weights in (2.4) all to be non- 
negative. One way to see this is to note that we can in this case decompose the node into 
a tree of 2 — * 2 nodes. The internal edges of this tree can be made to carry an arbitrary 
SU(2) matrix, which can however always be set equal to 1 by making a gauge transformation 
on the SU(2) matrices on the incoming and outgoing edges of the node (this is not always 
possible if there are cycles.) We may therefore introduce such matrices on each internal edge 
of the tree and integrate over them without changing the problem. Thus the weights for 
the node are products of weights for 2^2 nodes, which we have already argued are always 
non- negative. 

Next we consider whether this condition is necessary. Consider the terms in the expansion 
of det S which contain a factor S11S22 ■ ■ ■ Sn^^^n^^. The coefficient of this term is the 3x3 
minor det Sjj with I = J = {N — 2,N — 1, N}. The above Lemma about 3x3 matrices 
then shows that either this sub matrix has at least one vanishing element, or the product 
5*11 5*22 • • • 5'Ar_3,7v_3 Vanishes. In general, every 3x3 submatrix of S must have at least one 
vanishing element, or every term in the expansion of its complementary minor must vanish. 
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For A^ = 4, this implies that there must be at least 3 vanishing elements, not all in the 
same row or column. By considering the different cases, together with a suitably chosen 
Euler angle representation, it is possible to show that in each case a sufficient number of 
the 0(2) matrix elements must vanish that the diagram breaks up into a tree. This shows 
that, for A^ = 4, the condition of complete reducibility is also necessary for non-negative 
weights. However, we have not found a general argument for all A^ and indeed there may 
be exceptions. What can be shown straightforwardly is that S must have at least A^ — 1 
vanishing elements. 
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APPENDIX: SOME PROPERTIES OF SU(2) MATRICES 

We show explicitly that the integral 

1= I dU expizbiUbn) = 1 (Al) 

Any SU(2) matrix may be parametrised as t/ = exp{iaa ■ n) = cos a + ia -n sin a. The Haar 
measure is then 

I dU = (27r^)-i r sin^ ada f dn^ (A2) 

The exponent in (Al) has the form Acosa + in-B where A = zh\hR and B = zb\^abR. Note 
that B^ = A^. Although these are in general complex, since X is an analytic function of each 
of their components, we can first assume they are real. Then, without loss of generality, we 
can assume that B is real and points in the z-direction. Then 

X = (I/tt) / sin'^ ada (i(cos^) exp (A(cosq; + icos^sina)) (A3) 

Jo J-i 

The integral over cos^ is simple, and the result may be expanded in a power series in A. 
All terms, save that 0(A°), then vanish on integration over a. 

We show that any real linear combination of SU(2) matrices is itself, up to a real constant, 
an SU(2) matrix. From the above representation, it may be written as 

G = ^ flj cos ttj + ^ X! '^i^J ' ^ ^^'^ % (^^) 

i j 

which has the form A + iBN ■ a where A and B are real, and N is another unit 3-vector. 
Writing A = pcosa and B = psina then gives the required result. 
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